First order valence transition in Ce and YbInCu4 in the {B, T) - plane 
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The puzzling properties of the first order phase transition in YbInCu4 and its alloys in the wide 
range of magnetic fields and temperatures are perfectly described in terms of an entropy transition 
for free Yb ions. In particular, it turns out that the transition line in the (i3,T)-plane is very close 
to the elliptic shape, as it has been observed experimentally. Similar calculations are done, and the 
experiments are proposed for the (7 — a) phase transition in Ce in Megagauss fields. We speculate, 
that in case of YbInCu4 the first order transition is a Mott transition between a higher temperature 
phase in which localized moments are stabilized by the entropy terms in the free energy, and a 
band-like non-magnetic ground state of the /-electrons. 



PACS numbers: 75.30Kz, 75.30.Mb, 71.23.An 

For years the isostructural 7 — 0; transition in metallic 
Ce has been a classic example of the first order transition 
into a state with an intermediate valence (7 - Ce'' '^^+ 
a - Ce'^-^''+. For the phase diagram of Ce, seeEI. Change 
in the valence state was judged by the change in the 
unit cell volume). Although a structural transition in a 
crystalline matter possessing a critical point in the (P, T)- 
plane is of a great interest by itself, the discovery of a 
non-integer valence in a - Ce has opened the field of the 
so-called intermediate or mixed valence (MV) states in 
the rare earths and actinides (both for elemental metals 
and intermetallic compounds). 

The (7 — a) transition in Ce takes place in the pressure 
range P~10-^20kBar. Therefore, the isostructural tran- 
sition in YbInCu4H at r5'^40K and at ambient pressureS, 
looking akin to the Ce "isomorphic" transition, has at- 
tracted recently a lot of interest due to the possibility of 
study of that phenomena in much greater details. |-(Tihe 
major experimental results are best summarized inoB). 

In what follows we address the issue of the phase dia- 
gram of Ce and YbInCu4 in the {B, T)-plane, where B is 
a magnetic field, T is a temperature. Indeed, among 
many interesting results ofalj, the most suprising one 
is the universality of the first order transition line for 
YbInCu4 and its alloys. Namely, being expressed in the 
reduced variables {B/Bco,F/Tyo) the transition line sep- 
arating the high temperature phase (paramagnetic, local 
moments) and the low temperature "metallic" phase is a 
perfect circle (where T„o is the structural transition tem- 
perature in the absence of the magnetic field and Bco 
is the critical field at T=0).We will show that these re- 
sults are well described in terms of an entropy first order 
phase transition between the local /-moment phase and 
another phase probably of a less ordinary nature. An 
origin of this phase however seems not to be important 
if this second phase is characterized by a larger energy 
scale. The same ideas applied to the (7 — a) transition 
in Ce predict similar behaviour in high magnetic fields 
with Bco^200 Tesla. This is the field range achijevable 
for modern Megagauss magnetic field experimentgj. 



The valence of Cc in the 7-phase ^ very close to the in- 
teger /-occupancy (see, e.g. rewiewB), i.e., in the atomic 
configuration {Xe + 4/5(i6s^)t2l all d- and s-electrons of 
Ce go to the metallic bands. In accordance with the 
Hund's rule the ionic ground state has the total angu- 
lar momentum J = 5/2 which is split further in the 
cubic environment into a Fy doublet and a Fg quartet 
(in the Ce 7-phase F7 lies below Fg). Similarly, for Yb, 
the atomic configuration {Xe + 4/^'*6s^) results in the 
trivalent Yb^+ ionic configuration for the high tempera- 
ture YbInCu4-C15b phase, leading, to a localized /-hole. 
The f^^{J = 7/2) ground state is split by the crystal 
field into a quartet (F8)|-ajnd two doublets (Fg and F7). 
Inelastic neutron studiestll at T> 45K reveal the crystal 
field scheme with Fg and F7 lying at 3.2 meV and 3.8meV 
respectively above the ground state quartet Fg. 

The first order transition line is determined by: 

Fu{B,F) = Fl{B,T). (1) 

In (|l|) Fl and Fu stand for the free energies of the up- 
per and lower phases. The main assumption we use be- 
low is that the characteristic energies governing the be- 
haviour of the two phases differ significantly. We denote 
these scales as and T^, two effective "Kondo temper- 
atures" , in accordance with the existing tradition in the 
experimental literature fflifjlot data versus the isolated 
Kondo center propertieaaO'El (for extensive discussion of 
the theoretical results for the degenepate Anderson mod- 
els and the experimental results, seella). 

For YbInCu4 ~ 2bK while ~ 500Miy. With 
Ty , the temperature of the "valence transition" , for Yb 
and its alloys-Jying in the range of 10 — \QQK and 
SOTeslaaJa. The Fl{B,T) in (|l]) can be taken as 
a constant, neglecting the magnetic susceptibility term, 
while for the Fu{B,T) with the trivalent Yb'^^ consid- 
ered as a local free ion, one has: 

Fu{B,F) = Eo-T-S{B,T) (2) 

with the band energy Eq being actually temperature in- 
dependent below (this assumption is discussed in more 
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details below). Correspondingly, the first order transition 
line in the {B, T) plane is given by the equation: 



T • S{B, T) = const, 



(3) 



where the entropy is determined by the Yb^+ multiplet 
structure only. 

The magnetic susceptibility x{T) of YbInCu4 above 
r„o = 42i4r follows the Curie- Weiss law with an effective 
moment only negligible (by 5%X-smaller than the whole 
J = 7/2 ground state momentD'Ej. Thus, we first neglect 
the crystal splitting and write: 



T-S{B,T) = -Tin I ^ exp 
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where gj is a g-factor (for J — 7/2,gj ~ 8/7). From 
(^ the relation ansBco = T^o between the critical field 
Bco at T = and the value of the structural transition 
temperature T^o at zero field, is of the form: 



gjJfiBBcQ = 11,0 In (2J + 1) 



(5) 



which gives a « 1.9 for J = 7/2, a result whici is re- 
markably close to the experimental value a ~ 1.83. 
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FIG. 1. Function R{4>) (see text, Eqs.(9,10)) calculated for 
Yh^+ (J=7/2). Deviations of i?(0) from 1 do not exceed 6%. 



We re-write the Eqs. (|^,^ for the phase transition 
line using the new reduced variables /3 = B/Bco and 
T = T/TyQ and with the help of (ph we obtain: 




-m ( ;^ln(2J+ i; 



(6) 



Using the parametric form /3/r=tan0 and the identity: 
^2 + = cos-^ (j), (7) 
one may re- write (o) as: 



(8) 



where 

R{(j)) = In^ (2J+ 1) X 
J 



cos 6 ■ In 
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The plot of the function R{(f)) is shown in Fig.l. 

Since the deviation of R{(t)) — 1 from zism, does not 
exceed 0.06, we arrive to the main result ofEH13: 
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Postponing a detailed discussion for further publication, 
it is nevertheless necessary to highlight some essential 
features of our picture. The first interesting question 
is whether the account of crystal field split multiplets 
would improve the overall agreement with experiments. 
Although we have analyzed the relation: 



a^BBco — 71,0 
in terms of the crystal field Hamiltonian: 

H = H crystal + gjlJ-BJ ' B 



(11) 



(12) 



we do not stay on the results here, because (11) must 
display some cubic anisotropy which was not experimen- 
tally studied yet (two components in (^2|) do not com- 
mute with each other). We will limit ourselves with a 
comment that the energy levels' scheme for ( [l^ ) follows 
straightfo*wardly from making use of the explicit wave 
functionsEil for the representations r6,r7 and Fg. For 
the magnetic field, applied along the main cubic axis, it 
turns out that the experimental value a ~ 1.8EII1j is again 
closely reproduced in such analysis. We also would like 
to emphasize that the entropy: 



^(r,o)^ln<i4 + 2exp( 

J-vO 



2exp(-|Z-)^, (13) 



with Eq,E7 taken from0, is rather-, |elose to its value 
0.8 In 8 as integrated through T^o (seeQ'&y) for YbInCu4, 
which indirectly confirms the applicability of an isolated 
crystal field split hole state for Yb'^^ paramagnetic ion. 

The Yb"^+ hole occupation in the high temperature 
state determined from Yb-Ls X-ray absorbtion for the 
most compositions studied incl turns out to be really close 
to the Yb^+ trivalent state. This last fact, however, 
does not preclude yet that the upper phase may have 
developed pronounced Kondo, (effects with Tk — 25K, 
as e.g. is stated in referencesB'tHl^. On the other hand, 
it is not clear whether the existing data show consider- 
able deviations from the free ion behavior for the up- 
per phase. However, if it were so, (||) would not have 
been correct at low temperatures. In such case, one may 
choose for Fu{B,T) another expression, say, the exact 
solution for the exchange model or for the degenerate 
Anderson model. In the non-magnetic phase the scale, 
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is rather large and one may neglect the 
temperature dependence in Fl {B, T) at temperatures be- 
low 50— lOOiC. As for the conduction band electrons, typ- 
ical energy scale for Ce would be of order of leV. Such a 
scale for YbInCu4 and its alloys comprises probably only 
~ O.leV, as discussed below. 

It would be interesting to check, of course, whether the 
circular shape (|o|) of the transition line in the {B,T)- 
plane is indeed due to the entropy transition in the free 
ion scheme of Eqs.(^,^ or the result could be merely 
robust numerically. Unfortunately, the Anderson model 
thermodynamics in high magnetic fields has been studied 
in the Coqblin-Schriffer model limit (the charge is fixed) 
for Cc (J = 5/2) but not for Yb (J = 7/2). Even for 
Ce, there are published r*esults only for magnetization 
and specific heat (see inlla). To obtain the free energy 
expressions, one would need to integrate these data back, 
or solve the Bethe Anzats equations again. 
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FIG. 2. The line (solid) of the first order phase transition in 
the {B, r)-plane for Ce 7 — a transition. In the reduced vari- 
ables, deviations from the perfect circle (dashed line) would 
not exceed 6%. 

In Fig. 2 the phase diagram in the {B, T)-plane shows 
the first order phase transition line (solid), calculated 
for Ce according to (^ (J = 5/2, gj — 6/7). Its shape 
is again close to a cirlce in reduced variables (dashed 
line): the deviations from the circle do not exceed « 8%. 
The metamagnetic (7 — a) transition in Ce has not been 
measured yet, to the best of our knowledge. One sees that 
the low temperature values of magnetic fields are in the 
experimentally accessible Megagauss range {B"' ~ 185 
Tesla). The experiment presents a considerable interest 
and allows one to verify the free ion model we are using. 

Finally, let us discuss the physics, which may be re- 
sponsible for the transitions in Ce and YbInCu4. A first 
attempt to describe the Ce (7 — a) tcaiisition was the 
Falicov-Kimball-Ramirez (FKR) modelta. Although the 
FKR model is capable to reproduce the appearence of the 
critical point on the (7 — a) transition line in the (P, T)- 
plane, it does not reproduce such crucial feature of the 
a- phase as its intermediate valence. 

So far, the (7 — a) transition in Ce is usually inter- 



preted-in-|terms of the Kondo Volume Collapse (KVC) 
modetj'EZl. In the KVC model Ce atoms at the transition 
are treated as Ce'^+-ions in the both a and 7 phases (ap- 
proximately one electron in the /-shell), although in the 
two different Kondo regimes. As it is known, the Ander- 
son impurity model reproduces the Kondo behaviour in 
the regime when charge fluctuations are fully suppressed, 
and provides for the Tk the expression: 

^KOcexp|-^^|, (14) 

where | Ey- | is the effective position of the localized level 
below the chemical potential and the level's width F cx; 
V^v{tp) depends on the hybridization matrix element,!^, 
and the density of states at the Fermi level, v{eF)- 

The KVC modelta connects the first order transition 
with strong non-linear dependence of the Kondo scale 
(0) (I £y 1^ r) on the volume through the volume de- 
pendence of the hybridization matrix element (change in 
the volume of a unit cell is large, 5v/v ~ 20%!). Our 
arguments in the beginning of the paper, regarding the 
two different Kondo scales as needed for the applicability 
of (|^), agree well with the valuesEZi: 

T^~81.2± 12.2mey; r;^~8.2 ± 1.5mey 

Nevertheless, the KVC model seems not to be appli- 
cable in case YbIriCu4, where the volume changes are 
extremely smallatTLl. On that-xeason, the FKR model 
has recently been revisited inEEl. It is interesting, that 
although being somewhat sensitive to the choice of the 
model parameters, the elliptic shape of (11) for the phase 
transition Imp in the {B, r)-plane is preserved in the 
calculationgll3. This is probably due to the same mech- 
anism as above , i.e. due to large differencies in the en- 
ergy scales for the two phases (it seems however that 
the constant a in (12) strongly depends on the param- 
eters choice). Nevertheless, the FKR model can hardly 
be applicable for the YbInCu4 compound. In addition 
to its well known drawbacks, such as an absence of hy- 
bridization, large changes in the Uf occupancy, it seems 
that the peculiarities of this compound may originate in 
somewhat unusual features of its non-magnetic analog 
LuIiiCu4. This point has already been discussed inta. 
IrO the authors suggested a prechanism based on the 
band structure calculationsE2lEil for the semimetallic state 
observed both in LuInCu4 and YbInCu4 (Yb^+ , the high 
temperature phase!)Eli3. In this state Sf level falls into a 
"pseudogap" (or a dip in the DOS) at the choice of the 
chemical potential corresponding to the Yb'^^ configu- 
ration. If the exponential form for Tk in (|lj) remains 
correct, the strong non-linearity in (|l^ ) comes from the 
rapid—changes in the values of the DOS at the Fermi 
ZeweO. Thus, the high temperature state (the one above 
Tt,o — A2K) is stabilized by the entropy gain, while the 
phase with the higher DOS at the Fermi level is prefered 
at lower temperatures. The new state would, hence, coip 
respond to a non-integer valence Yb^'*+ as measured inEl 
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to be Yb^'^+ (an interaction term needs to be added to 
stabilize the valence). 

The feature, which remains not well understood with 
the above explanation is that the Hall coefficient sharply 
decreases in the low temperature state, as if the num- 
ber of carriers is getting comparable with the stoichio- 
metric value for the divalent Yb-ion (no hole in the f^'^- 
shell). The increase-|is too large compared with the va- 
lence change ~ 0.2EI. The energy scales involved and 
change in the [Hall coefficient are also not consistent with 
the results ofEJO. Indeed, AE, the energy change per 
Yb using Eqs. (|^,|) is: 

A£; = T„o In (2 J 1) - 90/^" (15) 

i. e. is too small to account for the large variation in 
the number of carriers. We propose another view on this 
problem, namely a weak Mott transition. At T > T^q 
localized moments Yb'^+ are stable due to the entropy 
gain, and exist as the localized holesli3. (Yb'^+)CuIn4 
is a band-like semimetal with a small carriers concen- 
tration. Below TyQ the valence change is small because, 
unlike in the FKR model, only interactions with small 
electron-hole pockets are essential. We speculate that af- 
ter the transition into the low temperature phase even 
the /-electrons form a band state, so that a small change 
in "occupation" numbers does not contradict to an emer- 
gence of a large /-like Fermi surface. An indirect support 
to these views mstf be fojtuad in the recent band struc- 
ture calculationsE2l. Fromc3, one may conclude that the 
change of the Yb valence ~ 0.2 at the 42K would just 
result in a shift of the chemical potential by ~ 0.01 eV 
inside the strongly featured DOS-sttith a broader width of 
the order of O.leV (see Fig. 8 inE3). This shift provides 
a correct magnitude for the energy change as estimated 
by (|l5|). The band picture is also in a reasonable agree- 
ment with rather high value for the Sommerfeld coeffi- 
cient in the linear term for the electronic specific heat, 
7 « hbmJ/mole ■ . 

To summarize, we have shown that the entropy transi- 
tion between the free ion paramagnetic state and the low 
temperature metallic state perfectly explains not only the 
elliptic shape of the transition line in the (S, r)-plane 
but also provides correct numerical results for its param- 
eters. Basing on similar calculations, we have suggested 
an experiment on the metamagnetic transition between 
the (7 — a) phases in cerrium. As for the true nature of 
the transition itself, we suggest that in YbInCu4 it is a 
weak Mott transition between a /-band metal and the 
semimetallic phase with the localized Yb'^+-holes. 
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